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Abstract 

Statistical and phase properties and number-phase uncertainty relations 
are systematically investigated for photon states associated with the Holstein- 
Primakoff realization of the SU(1,1) Lie algebra. Perelomov's SU(1,1) coher- 
ent states and the eigenstates of the SU(1,1) lowering generator (the Barut- 
Girardello states) are discussed. A recently developed formalism, based on 
the antinormal ordering of exponential phase operators, is used for study- 
ing phase properties and number-phase uncertainty relations. This study 
shows essential differences between properties of the Barut-Girardello states 
and the SU(1,1) coherent states. The philophase states, defined as states 
with simple phase-state representations, relate the quantum description of 
the optical phase to the properties of the SU(1,1) Lie group. A modified 
Holstein-Primakoff realization is derived, and eigenstates of the correspond- 
ing lowering generator are discussed. These states are shown to contract, in 
a proper limit, to the familiar Glauber coherent states. 
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I. INTRODUCTION 



The quantum mechanical calculations are based on appropriate sets of states in Hilbert 
space of a system. The basic system of quantum electrodynamics is the quantized single- 
mode electromagnetic field which is modeled by the quantum harmonic oscillator. The 
corresponding basic set of states is the complete orthonormal set of the number states |n) 
(n = 0, 1, . . . , oo), that can be used for expanding all the field states. From the other hand, 
the overcomplete set of the Glauber coherent states (CS) \a) has a number of remarkable 
properties and has been proved to be extremely useful in quantum optics The Glauber 
CS are closely associated with the boson creation and annihilation operators a) and a, which, 
together with the identity operator 1, form a realization of the Heisenberg-Weyl nilpotent Lie 
algebra J3J. The corresponding Lie group is the dynamical symmetry group of Hamiltonians 
for a number of important quantum mechanical problems [§,0. The Glauber CS |a) can 
be defined in three ways U[[7[: (a) the eigenstates of the lowering operator a; (b) minimum- 
uncertainty states or, more generally, intelligent states for position and momentum; (c) 
states constructed by action of displacement operators, which represent group elements, on 
the vacuum state. For the Heisenberg-Weyl group all these definitions are equivalent, but 
for other Lie groups they lead to distinct states. 

In the present work we concentrate on the SU(1,1) Lie group whose algebra has a number 
of realizations related to the quantized light field. The most known of them are the single- 
mode realization in terms of the amplitude-squared boson operators and the two-mode 
realization in terms of the boson creation and annihilation operators of the modes |5]||. 
These realizations and states associated with them have been studied in connection with the 
field quadrature squeezing [RMI1. We consider here the Holstein-Primakoff (HP) single- 
mode realization of the SU(1,1) Lie algebra [0,0. The analogous realization had been 
formerly introduced by Holstein and Primakoff [OJ for SU(2), and in the SU(1,1) case it is 
given by fl2| 



K + (k) = \]a)a + 2k-la\ K_(k) = a\]a)a + 2k-l, K 3 (k) = a^a + k. (1.1) 

Here k is the Bargmann index labeling unitary irreducible representations of the SU(1,1) 
Lie group [see text after Eq. ( |2.2| )]. Aharonov et al. [|T3[ have written the HP SU(1,1) 
realization in another form, by using the number operator h = a^a and the Susskind- 



Glogower exponential phase operators e 1 ^ and e~ % ^ Hl4| , |T5| . It was shown |[l^JT^J17|JlqJl^j20|1 
that the HP realization is related to the problem of the quantum description of optical 
phase. The generators ( |1 . 1|) of the HP realization with k = | have been used in the Jaynes- 



Cummings model Hamiltonians with intensity-dependent coupling [21]. The SU(1,1) Lie 



group is the dynamical symmetry group of these Hamiltonians. A multiboson version of the 



HP realization (|1 . 1| ) have been constructed by using generalized boson operators p2 |. 

Various states associated with the HP SU(1,1) realization exist in the harmonic oscillator 
Hilbert space. These states can be conveniently treated by using general group-theoretical 
techniques. One can consider the generalized CS obtained by action of SU(1,1) group el- 



ements on the vacuum state P^,p|,|6|, p!9|j20|] , and the eigenstates of the lowering generator 
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K_{k) (the so-called Barut-Girardello states | 2^ , |25| ). The present paper is devoted to the 
systematic investigation of properties of these states. It is known that the SU(1,1) CS are 
closely related to phase states [|T^JT6| , p~T|JT~8| , pr9|^0|1 . We proceed with a subsequent develop- 
ment of this relation by constructing a class of states characterized by simple phase-state 
representations [pj}[ . An especial attention is devoted to the study of phase properties and 
uncertainty relations between the number and phase observables. We start in Sec. |TJ with 
the SU(1,1) CS. Phase properties and number-phase uncertainty relations are examined 
using a recently developed formalism, based on the antinormal ordering of the Susskind- 
Glogower exponential phase operators [27,28]. In Sec. Ill we discuss in detail statistical and 
phase properties of the Barut-Girardello states |24],|2f|. These states have sub-Poissonian 
photon statistics and form an overcomplete basis in the harmonic oscillator Hilbert space. 
This basis was used to construct a diagonal representation of the density operator, which 
was shown to be well-behaved for nonclassical photon states [25|. A class of generalized 



philophase states (states with simple phase-state representations PS| ) is considered in Sec. 
[TV] , and states with sub-Poissonian statistical properties are found. A modified HP SU(1,1) 
realization is obtained by using the antinormal ordering of the exponential phase operators. 
Eigenstates of the modified lowering generator can be described as philophase states in one 
limit and contract to the Glauber CS in another. In Sec. [V] we summarize our conclusions. 



II. THE SU(1,1) COHERENT STATES 



The concept of CS has been generalized by Perelomov p3|fl|6[| and Gilmore f29|,|7| for 



arbitrary Lie group. Generalized CS are obtained by action of group elements on an extreme 
state of the group Hilbert space and thus can be created by using Hamiltonians for which 
given Lie group is the dynamical symmetry group. We start our discussion by a brief review 
of basic properties of SU(1,1) Lie group. For a given value of k, the operators ( |1.1| ) form a 
realization of the SU(1,1) Lie algebra [BU . 



[k 3 (k),k±(k)] = ±k±(k), [K4k),k + (k)\ = 2k 3 (k). (2.1) 

The group SU(1,1), whose generators are K±(k) and K 3 (k), is the most elementary non- 
compact non-Abelian simple Lie group. It has several series of unitary irreducible represen- 
tations: discrete, continuous and supplementary pOfl . In the present work we discuss only 
the case of the discrete series. The Casimir operator Q{k) for any irreducible representation 
is the identity times a number, 

Q(k) = [k 3 (k)] 2 - l -\k + {k)k_{k) + k_{k)k + {k)\ = k(k - i)i. (2.2) 

Thus a representation of SU(1,1) is determined by a single number k (the Bargmann index); 
for the discrete series this number acquires discrete values fc = ~,l,|,2,.... The correspond- 
ing state space is spanned by the complete orthonormal basis \k, n) (n = 0, 1, . . . , oo), 
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(k, n\k, n') = 5, 



nn' i 



\ k > n )( k > 



71 



(2-3) 



n=0 



These states may be defined by the following relations |3 

K 3 (k)\k,n) = (n + k)\k,n), 



K + (k)\k,n) = y /(n + l)(n + 2k )\k,n + l), 
K-(k)\k,n) = y/n(n + 2k - 1) \k,n - 1). 



(2.4) 



By considering the action of the operators K±(k) and K^{k) of the form ( |1 . 1|) on the number- 
state basis \n) of the harmonic oscillator Hilbert space, we see that this number-state basis 
coincide with the discrete series SU(1,1) orthonormal basis \k,n) for any allowed value of k. 
It means that the discrete series SU(1,1) state space is equivalent to the harmonic oscillator 
Hilbert space, when using the HP realization. Speaking elsewhere in the following about 
SU(1,1) we mean its HP realization defined by ( [Lip , with the orthonormal basis being 
merely the number-state basis. 



A. Coherent-state basis and the analytic representation in the unit disk 

As an example of the use of the group-theoretical methods, we consider the general 
results of Perelomov p3||5|,|6| for SU(1,1) discrete series CS. A coherent state is specified by 
pseudo-Euclidian unit vector t of the form 

r = (cosh r, sinh r cos tp, sinh r sin ip). (2-5) 

The CS \k,() are obtained by applying unitary operators £>(£) to the vacuum state, 

|*,C> = exp [tK + (k) - CK-W] |0> = (1 - |C| 2 ) fc exp [(K + (k)} |0>, (2.6) 

where £ = — \e~ %Lp and £ = — tanh^e - ^, so |£| < 1. Expanding the exponential and using 
~M, one gets the decomposition of the CS over the number-state basis, 



l^o = (i-ici 2 ) fc E 



n=0 L 



r(n + 2k) 
n\T(2k) 



1/2 



C\n). 



(2.7) 



The condition \(\ < 1 means that the 'phase space' of the SU(1,1) CS is the interior of the 
unit disk. The CS are normalized but not orthogonal to each other: 



i-KilT(i-|C 2 



2\k 



(i - crc 2 ) 2fc 

The identity resolution is an important property of the CS: 



(2.8) 



d»(k,()\k,()(k,t\ = i, 



(2.9) 
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where 



d(i(k,C) 



2k 
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;i-ici 2 ) 2 ' 



(2.10) 



and for k — % the limit k — > | must be taken after the integration is carried out in the 
general form. Thus the SU(1,1) CS form an overcomplete basis. 

One can can represent the harmonic oscillator Hilbert space as the Hilbert space of entire 
functions f{k,(), which are analytic in the unit disk. For a normalized photon state 



l/HEOX/Oln), 



n=0 



we get 



f(k,o = (i-\c\ 2 r k (k,e\f) = Y<Cn(f) 



n=0 



F(n + 2k) 
n\T(2k) 



1/2 



and this state can be represented in the coherent-state basis: 



|/>= / dfi(k, c)(i- \(\ 2 ) k f(k,C)\k, c>- 



(2.11) 



C, (2.12) 



(2.13) 



We will refer to such representations as the representations in the unit disk. The generators 
K±(k) and K^{k) act on the Hilbert space of entire functions f(k, () as first-order differential 
operators: 



d 



d 



d 



K +( k ) = ( 37 + 2fc C, K_{k) = -, K 3 (k) =(- + k. 



d( 



dC 



d( 



(2.14) 



B. Phase states and phase-like states 

We discuss now the CS \k, () with k = \. This case is interesting by two reasons. Firstly, 
phenomenological Jaynes-Cummings model Hamiltonians with intensity-dependent coupling 
have been constructed pi] ], f° r which the SU(1,1) group in the HP k = | representation is 
the dynamical symmetry group. In principle, the SU(1,1) CS with k = | can be created by 
applying such Hamiltonians to the vacuum state. Secondly, the SU(1,1) CS with k = \ are 
closely related to the phase states \9), given by |L5] 



1 oo 

|0> = ^=£^ n V>. (2-15) 



27T n=0 

For k = |, we can write, by using the number operator n and the exponential phase operators 

and e~ i( t>, 

K + {\) = ne 21 ^ = e^n, K 3 (i) =h+~. (2.16) 
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It was noted |16 , |13 ,17 that the CS 

00 

liO = v /l -Kl 2 Ec n |n) 

n=0 

are the eigenstates of the lowering exponential phase operator e 1 ^, 

e*liO = Cli,C>, 

just like the phase states, 

^|0) = e i ». 



(2.17) 



(2.18) 



(2-19) 



Therefore we will refer to the states ||,C) as the phase-like states. These states depend on 
two real parameters, |£| and arg£, while for characterization of phase only one real periodic 
parameter 9 is needed. When \(\ — > 1, one gets 



\0) 



hm(i-ici 2 r 1/2 iic). 



/2tt ICI-i 

Phase states \9) of Eq. Q2.15|) resolve the identity, 

rBo+2-K 



d6\6){0\ 



(2.20) 



(2.21) 



where 8q is a reference phase. Therefore, for arbitrary normalized state |/) of the form 
( |2.11| ), the phase-state representation is given by [^B] 



\f) 



'2n Je 



dee(f;9)\9), 



(2.22) 



where 



vnO 



Q(f;6) = V2^(6\f) = Y,C n (f)e-" 
The SU(1,1) generators act on the 0(/; 0) as first-order differential operators: 



e- w + ie~ ie 



d 

~dB' 



ka) = ie 
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d 



. d 1 



^2.23) 



(2.24) 



Thus the phase-state representation is a limiting case of the SU(1,1) coherent-state represen- 
tation when k = \ and |£| — > 1, that is, the ( representation is redefined on the unit circle 
instead inside. The function Q(f;8) of Eq. ( 2.23|) then can be called the "boundary func- 
tion" of the function /(|,C) analytic in the unit disk [Eq. ( |2.12| )1. The boundary function 
Q(f;6) determines uniquely the analytic representation /(|,C) p0| : 



n=0 



e ° +2 * de ®(f;e) 



00 



1 - (e ie 



(2.25) 
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In general, the Q(f;6) is given by Fourier series of the form (|2.23|) . However, there are 
states for which this Fourier series can be converted into a relatively simple functional form. 
We call such states the philophase states, and the phase-like states ||, () are an example of 
them: 



i-IC I 2 

1-Ce- 



9(i,C;0) = V2rc(e\l() = v -J-^ , (2.26) 



C. Statistical and phase properties 

The photon- number distribution of the CS \k,() is the negative binomial distribution 
H0: 



P n (k,C) = \(n\kX)\ 2 = (1 - iCl^^^^lCr". (2.27) 



The mean photon number (n)k,(; and the intensity correlation function g^l are [fLS 



ICI 2 

(^,C = 2A; T ^_, (2.28) 

(n} 2 fe)C " 1+ 2l- (2 - 29) 

The photon-number distribution P n (k,() is super-Poissonian > l), and for A; = | it 
becomes the thermal distribution |T^|T^]. 

Phase properties of a normalized photon state are obtained by calculating expectation 
values of appropriate phase-related operators. We define phase-related operators as operators 
that can be written in the form of Fourier-like series in the exponential phase operators 



27,28 



oo 

G = G 1 + E [ G n (e-^) + G- n (e«) } . (2.30) 



n=l 



Fourier coefficients are given by 

1 r 9o+2n . a 

G n = 7T d9G(9)e m9 } (2.31) 

27T J0 O 

where G(8) is a classical function corresponding to the operator G. It is easy to see that 



phase-related operators have diagonal phase-state representation [27 , 28 



G= / d9G(9)\9)(9\. (2.32) 
Je 
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An example of such operators is the Susskind-Glogower cosine phase operator C p^Jl5 
given by 

d= + [ 80+27T d9 cos 9\9)(9\. (2.33) 

2 v ' Je 

Functions G(9) must be 27r-periodic functions with convergent Fourier series. For a function 
G(9), which is not intrinsically 2n periodic, we must use the periodic expansion on the entire 
real axis, e.g., 

#per = 9 for 9 G [9q, 9q + 2tt) + 27r-periodic expansion on R. (2.34) 
Then Hermitian phase operator is given by [p7| , ffij[] 

4> = [ e ° +2n d9 9 pci \9}(9\ = e i + f; [e n (e^) n + L n [^) n ] , (2.35) 

J an 1 



1 f 6 ^ ... md ( ^o + vr, n = 0, 



where 

= oW d66elrM = \ T 'in9o " _z n ( 2 - 36 ) 

27T Je -^-e m0 °, n*0. 

u Km'' 

The expectation value of a phase-related operator G over a photon state is 

/■0O+27T 

(G) = Tr (pG) = / d9G(9)(9\p\9), (2.37) 

where p is the density operator of the state, and we have used the form ( |2.32|) for G. 

It is easy to see that such operators as C 2 , 2 , etc. cannot be written in the form ( |2.30| ) or 
( [2.32j ). In order to calculate the expectation values of these operators we use the antinormal 
ordering of the exponential phase operators e 1 ^ and e~ t( ^, defined as a procedure that places 
all raising operators e~ % ^ to the right of all lowering operators e 1 ^. The antinormal ordering 
restores the unitarity for the exponential phase operators, 

= l^e^l = 1, (2.38) 

where two * on either side of an expression are our notation of the antinormal ordering. By 
applying the antinormal ordering to the product of any two phase-related operators G and 
F, we obtain another phase-related operator 

IGFI= d9G(9)F(9)\9)(9\, (2.39) 

J 0o 

where both G and F can be written in the form ( [2.32j ) . Thus any two phase-related operators 
commute, 

l[G,F]: = 0, (2.40) 
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so, the phase variable is unique. The vacuum state may then be described as a state of 
a random phase, similar to all other number states. We have explained |28] that by the 
antinormal ordering we exclude nonrandom phase properties for the vacuum. From Eq. 
( ^■39j ), one gets 



lG n l= / d9[G(9)} n \9)(9\. (2.41] 

We obtain, for example, 



rd +2n 

:c 2 := de co S 2 e\e)(e\. (2.42) 



All the information about phase properties of a state is contained in the phase distribution 
function Q(9), given by 

Q{9) = (9\p\9). (2.43) 
Then the expectation value of a phase-related operator G is 

(G) = [ 90+27T d6G(6)Q(6), (2.44) 
Je 

and 

A A r 6 +2n 

CGFl)= d9G(9)F(9)Q(9). (2.45) 

It has been shown recently p8| that results for phase-related expectation values obtained 
by using the antinormal ordering are equivalent to those calculated in the frames of the 
Pegg-Barnett formalism |32| with proper limiting procedures [33]. For a pure normalized 
state \ f) of the form ( |2.11| ), the Q{6) function is given by 

1 1 00 

W;f) = |(e|/>r = =-|e(/;e)| a = — E C n (f)C* m (f)e- i(n - m)e - (2-46) 

Z7r Z7r n,m=0 

Often we meet a special case, in which 

C n (/) = |C , n (/)|exp(zn^ / ). (2.47) 

Then the phase distribution function has form 

1 00 1 f 00 1 

Q(f;9) = — E M lnl (f)exp[-in(9-^ f )}^—{l+2j2M n (f)cos[n(9-^ f )}\, 

27r n=-oo 27r I n=l J 

(2.48) 

where 

M n (f) = E |C m (/)C m+n (/)|, n > 0. (2.49) 

m=0 
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Substituting into Eq. ( |2.44|) this expression for Q(f; 0) and Fourier series representing G{6), 
we get 



(G) f = (f\G\f) = G + ^M n (f)[G n exp(-in? f ) + G„ n exp(in? f )}. 



(2.50) 



n=l 



Fourier coefficients G n are given by Eq. (|2.31|) . Expectation values are calculated for ob- 
servables represented by Hermitian operators. For a Hermitian phase-related operator G, 
the function G{6) is real, and G_ n = G* n . Then 



(G) f = G + 2 M n (f)Re [G„exp(-m^ / )]. 

n=l 

The expectation value of the Hermitian phase operator <p of Eq. ( |2.35| ) is 

OO 1 

($) f = 6 + vr + 2 £ -M n {f) sm[n(9 - <p f )\. 



(2.51) 



n=l 



The natural choice of the reference phase is 

% = iff ~ vr. 



(2.52) 



(2.53) 



Then the expectation value of the phase observable is at the center of the interval [9q, 8q+2tt): 



(<j>) f = iff = O +7T. 

The antinormally ordered square of the Hermitian phase operator is given by 



0O + 27T 



dO 6>p Cr 



where 



2tt Jbq 



d6 6 z e 



2 Jn9 



n=l 



§7T 2 + 27T0Q + 01, 71 = 0, 



(2.54) 



(2.55) 



in 



[2.56) 



With the choice fl2.53|) the expectation value of *0 2 * is 

7T 2 °° f— 1") n 



n=l 



Mn(/), 



(2.57) 



and the phase variance is 



vr 2 00 f-ll" 



n=l 



(2.58) 



In a similar way, we obtain 
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(C) f = Mi(f) cos ip f , 



(2.59) 



[l-M 2 (f)]+M 2 (f)cos 2 p f , 



(2.60) 



:(AC)j: = GO/ - (C)j = ^[1 - M 2 (f)\ + [M 2 (f) - M((f)} cos^ f . 



(2.61) 



When using the Susskind-Glogower sine phase operator S [|lj,y||, analogous formulas for 
(S)f, (lS 2 l)f and *(A5*)j* are obtained by replacing cos<^/ into sin^ f in Eqs. ( [2.59|) -( E2~6T 
Now we apply these general results to SU(1,1) CS \k,()- We get 

T(m + n + 2k)T(m + 2k) 1 1/2 



(I _ l^l 2 ^ 00 



r(2fc) 



m=0 



(m + n)\m\ 



ICI 



2;/( 



The simplest case is A; = |. Then .Mn(^,C) = |C| n > so ; 

l 1-KI 2 



27t1 + |C| 2 -2|C|cos(^-^ c )- 



Here 



(0)c = ¥>c = ar §C 



for all values of k. Also, we obtain 



:(A0) 



2 * 
2 ,l > 



7T 



+ 4E 



n=l 



n 



y-icr 



(C)i c = |C|cos^ c , 



:(ac)! 



^(i-ICI 2 ) 



(2.62) 

(2.63) 

(2.64) 

(2.65) 
(2.66) 
(2.67) 



The mean photon number (n)i^ = |C| 2 (1 — |C| 2 ) 1 shows that the quantum limit (small 
excitations) is obtained for |£| <C 1 and the classical limit (large excitations) corresponds 
to |C| -> 1. In the quantum limit, when |C| -> 0, *(A0)| * « vr 2 /3 and VAC)? * « |, 
that corresponds to the uniform phase distribution, i.e., Ci $) ~ (2tt)~ 1 . In the classical 
limit, when |C| -> 1, both *(A</>)| * and *(AC)| r * tend to zero, that corresponds to a 



state with a perfectly defined phase. In this limit the Q{\-, C; function is very narrow. By 
using the expression (|2.62|) for A4 n (k, (), we can calculate numerically the phase distribution 
function Q(k,£;0), given generally by Eq. (p. 48 ). The larger values of k are, the narrower 
the Q(k, (; 9) function is, and the better the phase of the state is defined, for given \(\. 
We can calculate numerically also the *(A0)|^* and *(AC)^*, by using Eqs. ( |2.58|) and 
( 2.61 ), respectively, for different values of k. The phase variance *(A0)|^* is independent of 
(pt = arg£. When |£| — > 0, the phase variance tends to the random value vr 2 /3, and when 
|C| —* 1, it tends to zero, for all values of k. The larger values of k are, the smaller the phase 
variance is, for given \C,\. A similar situation is with the cosine variance *(AC)^*. When 



ICI - 

zero. 



0, the cosine variance tends to the random value |, and when |£| 



1, it tends to 
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D. Number-phase uncertainty relations 



The number and phase operators form a Heisenberg pair of canonically conjugate ob- 
servables |34j . The number-phase commutation relation is given by p7| , |26| 



n. 



i(l-2n\9 )(9 \). 



(2.68) 



The additional term — 27rz|# )(#o| i n this relation takes into account the periodicity of the 
phase. Equation ( |2.68|) is a special case of the general commutation relation between any 
phase-related operator G of the form ( p.32j ) and the number operator n [|27l j26[l: 



[n, G] = i 



(10 

For the phase operator (f> of Eq. (|2.35| ), we get 

d 



(2.69) 



de 



9 per = 1 - 2n5(9 - 6 ) 



(2.70) 



and the relation ( 2.68Q follows immediately from Eq. (|2.69 ). The number-phase uncertainty 
relation reads 



(An) 2 :(A0) 2 :>l[l-27rQ(£ o )] 2 . 



We define the following function: 



(An) 2 :(A0) 2 ; i 

[l-27rQ(^)] 2 " 4' 



(2.71) 



(2.72) 



By evaluating this function, we can investigate the number-phase uncertainty relation (p.71[) 
for various photon states. 

The number variance for the SU(1,1) CS can be calculated from the function P n (k, () of 
Eq. ( ^27p . One gets 



(An)l , = 2A- 



ICP 



'*,c -(i-iciT 

By using the standard choice ( |2.53| ) for 8q, we obtain from Eq. ( |2.48 



^2.73) 



Q(k,(;6 ) = — 

Z7l 



n=l 



(2.74) 



where coefficients A4 n (k, () are given by Eq. ( 2.62Q . The phase variance can be calculated 
using Eq. (p.58|) . The simplest case is k — |. Then *(A0) 2 ^* is given by Eq. (|2.65|) , and 



Q(h Ci ^o) 



i i-ICI 

27Tl + |C| 



(2.75) 
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Then Eq. (gjg ) reads 



VG,C) 



1 



4(1 -ICI) 2 



2 oo (_i)n 

n=l ri 



^2.76) 



For \(\ — > 0, this function tends to the random-phase value vr 2 /12, while for \(\ close to 1, it 
blows up. Numerical calculations show that this limiting behavior is universal for all values 
of k. We see that, for small values of k, the CS \k, £) are far from satisfying an equality in 
the number-phase uncertainty relation ( [2.711) . However, for a large k, the function V(k, £) is, 
in an intermediate range of |£| values, close to its limit -r. The fact that the V(k, C) function 
blows up for | C| — > 1 is, at first sight, somewhat strange and contrary to the standard 
conception of "the classical limit" where uncertainties can be neglected. Indeed, the phase 
of the SU(1,1) CS becomes perfectly defined in the limit |£| — > 1. But the situation with 
the photon statistics is absolutely different. The relative photon-number uncertainty for the 
\k, Q states is 

An fe , c 1 



(n) k ,c V2k\C\' 

For \(\ — > 1, we obtain nothing similar to the classical zero uncertainty (unless k — > oo). 
For example, for k — |, the relative photon-number uncertainty tends to unity in the limit 
\C\ — > 1. This follows from the fact that the ||,C) states have thermal photon-number 
distribution, though they are pure states. 

One can also discuss the uncertainty relations for the number and cosine or sine phase 
operators. From the commutation relations 

[n,C] = -iS, [n,S}=iC, (2.78) 

one can deduce the uncertainty relations 

(An) 2 :(AC) 2 : > i(5) 2 , (2.79a) 



(An) 2 :(ASf: > ^(C) 2 . (2.79b) 
In order to investigate these uncertainty relations, it is convenient to define the functions 

s MM > i (2 . 80a) 

(S) 2 "4' v ' 



(An) 2 ;(Ag) 2 ; > 1 

2 ^ 2 - 4' 



(2.80b) 



For a pure state |/), satisfying the condition ( p.47| ), we get 
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M{{f )sm <p f 



and the 1Z 2 (f) is given by interchange of cosy?/ and sin^j. For iff = tt/4, the functions 
72-1 (/) and 7^ 2 (/) coincide and come to the function 



M(/) s (A n) r (AC) ? : + :(AS) ? : = (An)? 1 -^. (2.82) 



We return now to the CS 0. When = |, we get 

^ 1 (I,C) = [2(l-|C| 2 )sin 2 v , c ]- 1 . (2.83) 

For \(\ — > 0, this function tends to the value [2sin 2 <^] -1 > |, while for \(\ close to 1, it blows 
up. The SU(1,1) CS with other values of k behave similarly. When k is small, the CS \k, Q 
give a strong inequality in both uncertainty relations ( |2.79a| ) and (|2.79b[) . However, for a 



large k, the function U(k, () is, in an intermediate range of |£| values, close to its limit ~. It 
means that the SU(1,1) CS with large values of k approach to satisfy, for some |£| values, an 
equality in the number-cosine and number-sine uncertainty relations. The functions V(k, C) 
and U(k, () behave very similarly. Therefore, the number-phase uncertainty relations can 
be studied using the Hermitian phase operator (f> as well as using the cosine or sine phase 
operators. 

III. THE BARUT-GIRARDELLO STATES 



In this section we use results of Barut and Girardello (BG), who have constructed [24 
the eigenstates of the lowering generator K_(k), 

K_[k)\k,z) = z\k,z), (3.1) 

where z is an arbitrary complex number. The normalized BG states can be decomposed 
over the number-state basis, 

gk— 1/2 oo 

%z) = —j= V . In), (3.2) 

/l 2k ^(2\z\)n=oJn\T(n + 2k) 



where l v is the z/-order modified Bessel function of the first kind. The BG states are not 
orthogonal to each other, 

<Mi|M 2 > = ^-i( 2 v^) . (3 . 3) 
//2fc-i(2bi|)/ afc -i(2|^|) 



The BG states have a simple representation in the SU(1,1) coherent-state basis: 

Z(k,z,() = (1 - \C\ 2 r k (k,C\k,z) = = = exp«). (3.4) 

N // 2fc -i(2|«|)r(2A;) 

In the following discussion we will consider various properties of these states. 
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A. The identity resolution and the analytic representation on the complex plane 



In order to prove that the BG states resolve the identity, one must find the measure 
dfi(k, z) such that 



dfi(k, z)\k, z)(k, z\ = 1. (3.5) 

Writing dfi(k, z) = fi(k, \z\)d 2 z and integrating over the whole z plane, we find 

~ \n)(n\ [™ „ „ . IS \ z \ 2n+2k . . x 

27r E „JJ ■ u / d\z\»(k,\z\) T l 1 . . =1. (3.6) 



-n\T(n + 2k)Jo 1 1 ^ v " " hk-i(2\z\) 
By using the following formula p5| 



J™ dt t a K b (t) = 2^ 1 r (^P^j r ( 2 "4 ±i ) ■ Re (a ± 6) > -1, (3.7) 



where ^ is the i/-order modified Bessel function of the second kind, we see that Eqs. ( |3.5| ) 
and ( |3.6| ) are valid for 

dfi(k,z) = -K 2k -i(2\z\)I 2k - 1 (2\z\)<Pz. (3.8) 

7T 

Thus the BG states form, for each allowed value of k, an overcomplete basis in the harmonic 
oscillator Hilbert space. 

The harmonic oscillator Hilbert space can be represented as the Hilbert space of entire 
functions f(k,z), which are analytic over the whole z plane. For a normalized state |/) of 
the form (|2.11| ), we get 

z 1 y/n\T(n + 2k) 

and this state can be represented in the BG basis: 

/ 1/2 

\f) = [ d^(k,z) } Z} f(k,z*)\k,z). (3.10) 



The SU(1,1) generators K±(k) and K$(k) act on the Hilbert space of entire functions f(k, z) 
as linear operators pi[: 



7 72 J 

K+{k)=z, K_(k) = 2k— + z—, K 3 (k) = z— + k. (3.11) 

dz dz 2 dz 
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B. Statistical and phase properties and the number-phase uncertainty relations 



The photon- number distribution of the BG states \k,z), 

I r |2n+2fe-l 

P n (k,z) = \(n\k,z)\' 



I 2k -x{2\z\)nT{n + 2k) 



(3.12) 



is sub-Poissonian ||25|| . The expectation values for the number operator and its square can 
be easily calculated, 



hk(2\z\ 



{n)k, 



l 2k-l 



(2N) 



(n 2 )k,z = \z\ 2 -{2k- l)\z\ 



hk(2\z\) 
(2M)' 



'2fc-l 



The intensity correlation function is given by 

(n 2 )k,z- (n)k,z hk-i(2\z\) 



(2) 
9k,z 



<2fc 



(n)l >z h k {2\z\) 
In the quantum limit \z\ <C 1, we get approximately 



hk(2\z\ 



2k 

\z\ 



(2) 

9k,z 



2k 



2k + 1 



(3.13) 
(3.14) 

(3.15) 

(3.16) 



We see that the maximal antibunching g kz = | is achieved for k = |, while for large values 

(2) 

of k, the g k z approaches unity. Another interesting feature is that the intensity correlation 
function g kz of Eq. ( [3.16 ) is equal to the reciprocal of the coherent-state intensity correlation 
function g k ^, given by Eq. (|2.29|) , 



(2) 

9k,z 



(2) 

9k,c 



-i 



\z\ < 1. 



(3.17) 



(2) 

It can be verified that the g k z is always less than unity, so the BG states have sub-Poissonian 
photon statistics. Therefore, we call these states (which are, in many aspects, similar to 
the Glauber coherent states) the subcoherent states. The detail discussion of statistical 
properties of the BG states is given in Ref. p5 |. Also, it is shown there that the BG basis 
can be used to construct a diagonal representation of the density operator (the so-called 
subcoherent P-representation), which is shown to be well-behaved for nonclassical photon 
states. 

According to the general results of Sec. [II C| , we investigate here phase properties of the 
BG subcoherent states. By using Eq. (|2.49|) , we find 



M n (k,z) 



\n+2k— 1 oo 



1 2m 



E 



I 2k -i(2\z\) ^ [m\(m + n)\T(m + 2k)T(m + n + 2k)} 1 / 2 ' 



(3.18) 



Then the phase distribution function Q(k,z;6) can be calculated from Eq. ( p.48j ). In the 
simplest case k — |, we get 
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and 



Here 



I\ n \(2\z\) exp[-in(0 - (p x )]. 



27T/ (2|Z|) „ 

ip z = argz for all values of k. By using the expansion [35 

oo 

e xcosd = £ I|n|(x)e- me , 



(3.19) 



(3.20) 



(3.21) 



we can write 



Q(f,z;0) 



exp[2|2;| cos(0 — <p z )] 



27r/ (2|z|) 

This result can be also obtained by noting that 



0(1 z;9) 



Io(2\z\) n=o nl y/lo(2\z\) ' 



(3.22) 



(3.23) 



We see that the BG subcoherent states \^,z) are an example of philophase states. In 
the next section we will show that the states \\,z) are a special case of a wide class of 
philophase states, all of them have the phase distribution function of the form fl3.22j ). This 
phase distribution function is the same as the classical von Mises distribution |]36fl . In the 
quantum limit |z| <C 1, the Q(^,z;8) tends to the uniform phase distribution, while in the 
classical limit |z| 3> 1, it is narrow. Also, we calculate 



:(A0) 



2 * 



7T 



+ 



11 y 

h(2\z\) 
W\A) 



-l) n 



In{2\z\ 



COS (p z 



(3.24) 



(3.25) 



:(AC)f 



2 * 

z* 



1 - 



h 


[2 


A) 


h 


[2 


A) 



+ 



h{2\z\) 


il(2\A) 


h{2\z\) 





In the quantum limit \z\ <C 1, we get approximately 

„2 



:(A0) 
:(ac) 



7T 



4UI 





M 2 







cos 2 <£> 2 



cos (p z . 



(3.26) 



(3.27) 
(3.28) 
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Thus the phase variance *(A0)f * and the cosine variance *(AC)\ * tend in the limit 

2 ' z 2' Z 

\z\ — > to their random values 7r 2 /3 and |, respectively. In the classical limit of large 
\z\, the *(A0)i * tends to zero, and the *(AC)f * behaves similarly tending to the value 

2 ,2 2 ,Z 

(2|z|) _1 sm 2 {p z . Thus, for large values of \z\, the phase of the states ||, z) is well defined. By 



using the expression (|3.18|) for A4 n (k, z), we can calculate numerically the phase distribution 
function Q(k, z, 8) for different values of k. The larger values of k are, the flatter the Q(k, z, 8) 
is, and the worse the phase of the state is defined, for given \z\. This behavior is opposite to 
that of the SU(1,1) CS. We calculate numerically also the variances *(A0)| 2 * and *(AC)| 2 * 
for different values of k. The phase variance *(A0)f, 2 * is independent of (p z = axgz. When 
\z\ — > 0, the phase variance tends to the random value 7r 2 /3, and when \z\ is large (\z\ ^> k), 
it tends to zero, for all values of k. The larger values of k are, the larger the phase variance is, 
for given \z\. We meet a similar situation for the cosine variance *(AC)| 2 *. When \z\ — > 0, 
the cosine variance tends to the random value |, and for large \z\ it tends to zero. 

Now we discuss the number-phase uncertainty relations for the BG subcoherent states. 
The number variance (An) 2 . z can be deduced from Eqs. (|3.13|) and ( [3.141 ), 



(An)L=M 2 -(2*-l)|*| 



hk(2\z\) 



|2 Il(2\z\ 



'i»_i(2|z|) 1 1 / 2 2 fe _ 1 (2|z|) 
By using the standard choice (|2.53|) for 8 , we obtain from Eq. Q2.48Q 



(3.29) 



Q(k, z; 8 



2tt 



l + 2j2(-l) n M n (k,z) 



n=l 



[3.30) 



where the coefficients Ai n (k, z) are given by Eq. ( |3.18|) . The phase variance can be calculated 
using Eq. ( p.58j ). The simplest case is k = \. Then *(A0)| 2 * is given by Eq. ( |3.24j ), and 



Q(k,z;8 ) 



exp(— 2\z\) 



27t/ (2|z|) • 

Thus we find the V function, defined by Eq. Q2.72f ), 



VQ,z) = \z\ 



1 ggN) 



Jn 2 (2kl 



1 - 



exp(— 2\z\) 
2uh{2\z\) 



n 2 A 00 (— 1) n 



3 I (2\z\) ^ 



rr 



In the quantum limit \z\ 1, the function V(|, z) is given approximately by 



(3.31) 



(3.32) 



V(|,z)« 



7T 
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Ul fi 



7T" 



12 



(3.33) 



In the classical limit \z\ ^> 1, the V{\,z) tends to its minimal possible value \. Numerical 
calculations show that this behaviour is universal for all values of k. The larger values of k 
are, the slower the function V(k, z) approaches, as \z\ increases, its limit 4. 

An analogous situation is with the number-cosine and number-sine uncertainty relations. 
For k — we get 
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2 ' 



sin 2 y9 2 



mm) 



1 - 



W\z\) 
W\z\) 



+ cos tp z 



W\z\) I*(2\z\) 
P (2\z\) 



(3.34) 



and the TZ 2 (h,z) is obtained by interchanging cosy^ and sin<y2 2 . In the quantum limit 
\z\ <C 1, we find 



n x {\,z) « (2sinV 2 )- 1 > - 



(3.35) 



while in the classical limit \z\ ^> 1, the result is 



(3.36) 



We see that, for large \z\, the 7£i( 



2 ' 



approaches its minimal possible value \. For other 



values of k we meet a similar behaviour. The larger values of k are, the slower an equality 
is achieved, as \z\ increases, in the number-cosine and number-sine uncertainty relations. 
Generally, the number-phase uncertainty relations for the BG subcoherent states behave 
very differently from those ones for the SU(1,1) CS. However, we again see that the functions 
V and U behave very similarly. Hence, it is not very important what is the phase function 
that we choose for investigating number-phase uncertainty relations. This fact is a direct 
result of the unique phase definition in the antinormal ordering formalism. We can conclude 
by noting that the subcoherent states \k,z) (especially ||,-z)) are, for \z\ 3> k, a good 
approximation to the number-phase intelligent states, i.e., they tend to satisfy an equality 
in the number-phase uncertainty relation. This property remains valid also regarding to the 
number-cosine and number-sine uncertainty relations. 



IV. PHILOPHASE STATES 

In the preceding sections we have seen that the states 1 1, an d 1 1, z) are two examples of 
philophase states, i.e., states for which the phase-state representation function 0(6*), defined 
by Eq. fl2.23| ) as Fourier series, can be converted into a relatively simple functional form. 
The function 6(§,C;#) of the SU(1,1) CS ||,C) is given by Eq. (£2§), and the function 
Q(^,z;9) of the BG subcoherent states \\,z) is given by Eq. ( |3.23| ). The philophase states 
[phase-like SU(1,1) CS] are the eigenstates of the exponential phase operator e 1 ^ [see 
Eq. ( |2.18| )1, and the philophase states \\,z) (BG subcoherent states) are, by the definition 
( p.l|) , the eigenstates of the lowering generator K-{\) = e^n. We can consider a generalized 
operator, defined by 

Z{a) = AT_(i) + ae^ = e^(n + a), (4.1) 

where a is, generally, any complex number. The eigenstates of this operator are obtained 
from the following eigenvalue equation: 

Z(a)\z,a) = z\z,a), (4.2) 
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where z is an arbitrary complex number. The states \z,a) can be decomposed over the 
number-state basis: 

oo 

\z,a) = Y,Cn(z,a)\n). (4.3) 

n=0 

From the eigenvalue equation ( |4.2j ) we deduce the following recursion relation for the 

C n (z,a): 

C n+1 (z, a) = Z C n (z, a), n>0. (4.4) 
n + a + 1 

The solution of this equation is 

z n 

C n (z, a) = C (z, °) T ( n + (T + 1 y (4-5) 

where Cq (z, <j) remains to be determined from the normalization condition. In the following 
discussion we restrict ourselves to the relatively simple but important case of integer o. 
When a is zero, we obtain the BG subcoherent states whose properties have been 

discussed in detail in the preceding section. It is convenient to distinguish between positive 
and negative values of a, since properties of the \z, a) states are essentially different in these 
two cases. 



A. The case of integer a < 

For a nonpositive integer a, the normalized coefficients C n (z,a) are 

1 z n-\er\ 

n n > W\ 



CSr ) {*,*) = \ ^o(2N)("-M) ! ' ' ' (4.6) 

0, n < \cr\. 

The index ' — ' stands here and in the following for negative values of a, though all results 
are valid also for a = 0. It is interesting to note that the states 

i 00 7 n -W\ 
M- = I , ■ „ E 7 i-ml n > ( 4 - 7 ) 



are not only the eigenstates of the Z{p) with eigenvalues z, but also the eigenstates of the 
operator 

t{z) =n- ze-** (4.8) 

with eigenvalues \a\ = —a, 

£{z)\z,a)_ = -a\z,a)_. (4.9) 
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We deduce from this equation the following recursion relation: 



n + a + 1 



Ct\z,a), n>-l, C { _- 1 \z,a)=0. 



(4.10) 



The difference between recursion relations flll]) and Q47T0D is very essential, since Eq. ( |4.4|) 
is meaningful for arbitrary values of a, while Eq. (|4.10| ) has a nonzero solution only for a 
nonpositive integer a. 

Equations ( |4.9| ) and ( |4 . 1 0| ) can be transformed into the differential equation [ f26|j : 



'. d 

l d9 + a ~ Ze ~ 



0. 



(4.11) 



Here, by the usual notation, >(z, a; 9) is the phase representation function for the states 
\z,a)-: 



0H( 



^(9\z,a)- = J2 C t\z,<7)e- in9 , 



(4.12) 



n=0 



\z,a)- 



1 



'2-k Je 



9 +2n 



d90 ( -\z,a;6)\9). 



(4.13) 



The normalized solution of Eq. (|4.11| ) is 



Q ( -\z,a;9) 



exp (ia9 + ze ld 



I (2\z\ 



(4.14) 



The Q('(z, a; 9) function must be 2tt periodic. Therefore a must be an integer. Moreover, 
from Eq. ( |4.12| ) we see that the Q(~\z, a; 9) function written as Fourier series must have 
only Fourier coefficients with non-negative values of n (this is a general requirement to the 
phase representation functions p6[). For the function ( |4.14|) this demand forbids positive 



values of a. From Eq. ( 4.14 ) we see that the states \z, er)_ are philophase states, and the BG 
subcoherent states \\,z) are a special case corresponding to a = 0. The \z, cr)_ states with 
different values of a have the same phase properties because the phase distribution function 



Q 



Hi 



2tt 



|e<->(; 



exp[2|z| cos( 



2nl (2\z\) 



arg z 



(4.15) 



does not depend on a. This function was considered in Sec. [Tl [B[ There we have discussed 
in detail the phase properties of the states \k = ~,z) = \z,a = 0), and this discussion 
remains unchanged with regard to all the states \z,a)~. 

The statistical properties of the philophase states \z,a)- can be calculated by using the 
photon-number distribution 



P<r\z,a) = \(n\z,aU 





\z 


2(n-H) 




*\) 




tr\)1\* 



n > lerl 



(4.16) 



n < \a\. 
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Another way is associated with the phase-state representation. For a pure state |/) with 
the phase representation function Q(f;6), the number-operator moments are given by [p6j 



In either way, the calculation is simple, and we obtain 

^i(2M) 



(n) 
(n 2 ) 



(-) 

z.a 



a. 



2\H 
z. a 



a 2 + \zr - 2o\z 



h(2\z\) 
h(2\z\) 



(4.17) 

(4.18) 
(4.19) 



The number variance, 



(An)^(-) = |*| 



If 


2 




I 2 o 


;2 





(4.20) 



is independent of a, just as the Q(~>(z; 0) function, and therefore the number-phase uncer- 
tainty relations for the philophase states \z, a)- are the same as those ones for the BG states 
\k = ~, z) = \z, a — 0). We have discussed these number-phase uncertainty relations in Sec. 
[JIB. The intensity correlation function is given by 



gfli 



<J + \z\ 



(2(7 + 1)1* 



/o 2 (2M) 



a' 



2a\z\ 



I?(2\z\) , , , 2 J?(2|z|) 



+ z 



P (2\z\) 1 1 J 2 (2|*|) 



-i 



(4.21) 



The gfl(— ) is always less than unity, so the photon number distribution \z,a) of the 
philophase states \z, er)_ is sub-Poissonian. For \z\ <C 1, we get 



a. 



(n 2 } { z ^ - (r 2 + (I - 2a)\z\ 2 + a\z\ 
9$U - 1 





* 


4 /2 




a 






a 2 - 2(7(1 


z 


2 _ 


* 


72) + 


* 


4 



(4.22) 
(4.23) 

(4.24) 



For a = 0, we obtain g Z Q « |, in accordance with Eq. (|3.16| ). When |a| > 1, expression 
( [4.24| ) can be further approximated: 



&{-) 



1 - 



CM 



(4.25) 



This result is interesting since, for a = — 1, the intensity correlation function tends to its 
minimal allowed value: g^ 2 \(— ) ~ 0. This is the maximal accessible antibunching. For 
a = —2, the gi 2 l(— ) again approaches |, as for cr = 0. With further increase of the 
(yfp^(— ) tends to the Poissonian value 1. The intensity correlation function ( |4.25| ) is the same 
as that for the number states \n) with n — \a\. It is not surprising, because we see from 
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Eq. 



that the philophase states \z,a)- tend in the limit \z\ — > to the number states 



\n 



\a\). For \z\ ^> 1, we get 



(-) 



(h) 



1^1 — °" 
(X 2 + |z| 

1 



i 
i- 



2(j(|2| 



1), 





z\/2 


<x 2 + 


z 


2 -2a{\z 


-i)- 


*|/2 



When 1 2 1 ^> |er|, expression (|4.28 ) can be further approximated: 

1 



1 



2\z\ 



(4.26) 
(4.27) 

(4.28) 



(4.29) 



All the BG states in the classical limit behave according to Eq. ( f4.29|) , and this behavior 
also holds for the eigenstates of the generalized operator Z(o). 

We finish the discussion of the case a < by noting that the philophase states \z,o)- 
arise from the problem of finding intelligent states for operators which are combinations of 
the number and phase- related operators ||26|| . We define two Hermitian operators: 

(4.30) 
(4.31) 



t x {z) = n~ (Rez)C - (lmz)S, 
± 2 {z) = (Rez)S- (lmz)C, 



such that 



£(z) = ±i(z) +i± 2 (z). 



(4.32) 



Therefore, the eigenvalue equation ( O ) means that the states \z,a)- are intelligent states 
for operators £1(2) and £2(2), i-e., for these states the uncertainty relation 



[:(As 1 ) 2 :][:(As 2 ) 2 :]> i ic[Si,E 2 ] 



(4.33) 



is an equality. Indeed, we can use results of this and preceding sections [Eqs. ( |3.25| ), ( |3.26|) 
and ( g7Tg )— ( gTgg )] in order to calculate 



z\h(2\z\) 



The commutation relation is 

l[t u t 2 ]: = z[(Rez)C + (lmz)S]. 

The expectation value is easily obtained by using Eq. (|3.25|) : 

h(2\z\) 



2 / (2|z|)- 



GPi,s a ]:)i-> 



i\z\ 



(4.34) 



(4.35) 



(4.36) 



and it is evident now that the states \z,a)- are the intelligent states for operators Si(-z) 
and i^U). 
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B. The case of integer a > 

As soon as a passes through zero, the properties of the states \z, a) sharply change. For 
a > these states are 

M+ = 1= f, ■ ; E 7^^»' ( 4 - 37 ) 



T (\z\,a) n=o 



[n + a ! 



where we have defined 



(7—1 



To(kU)^I (2N)-E£n2- ( 4 - 38 ) 



m =o l™ ! ) 2 ' 



The index '+' stands here and in the following for positive values of a, though all results 
are valid also for a = 0. The states \z, a) + are normalized but not orthogonal to each other: 



+ (zi, a\z 2 , a)+ = , ■ (4.39) 

IT (\z 1 \,a)T (\z 2 \,a) 



In the following we will consider a modified HP realization of the SU(1,1) Lie algebra. After 
that we will discuss statistical and phase properties of the \z,a) + states and their relation 
to the Glauber CS. 



1. Modified HP realization 
It is easy to prove that the states \z, a) + resolve the identity: 

J dfx(z,a)\z,a) ++ (z,a\ = 1, (4.40) 

where 

dfx(z,a) = -K (2\z\)T (\z\,a)d 2 z. (4.41) 

7T 

Indeed, Eq. (|4.4q) reads 

4E rJ W | ( 1,12 rd\z\K (2\z\)\zr^ = i. (4.42) 

By using formula (p.7|) , we make sure of this equality. Thus the states \z,a) + form, for 
each integer a > 0, an overcomplete basis in the harmonic oscillator Hilbert space. This 
result leads us to idea that the operator Z(a), whose eigenstates the \z, a)+ states are, is the 
lowering generator belonging to a realization of the SU(1,1) Lie algebra. This operator is a 
generalization of the lowering generator K^ik — |) = Z(p = 0), and therefore a modified 
realization with a > should be a generalization of the HP realization with k = |. We 
introduce the following sets of operators: 
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K_(a) = e**(n + tr) = Z(a), 

K+(a) = (n + a)<P* = [K_{&)}\ (4.43) 
K 3 (a) = h + a + \. 

These operators obey the SU(1,1) Lie algebra provided the antinormal ordering is used: 

:[K4a),K + (a)}: = 2K 3 (a), [K 3 (a), K ± (a)} = ±K ± (a). (4.44) 

[The subscripts '— ' and '+' of the SU(1,1) generators are not related to the sign of a] The 
antinormal ordering should be applied also to the calculation of the Casimir operator: 

Q = [K 3 {a)f - :\[K + {a)k_{a) + K_{a)k + {a)]l = (4.45) 

By comparing this result with Eq. (|2.2j) , we see that the modified HP realization Q4.43 ) 
corresponds to the case of the discrete series representation with k — \. The action of the 
generators ( |4.43| ) on the number states is given by 

K 3 {a)\n) = (n + a + %)\n), 

K + {a)\n) = (n + a + l)\n + l), (4.46) 
kS(a)\n) = (n + a)\n- 1). 



It follows from the comparison of these formulas with relations (|2.4|) that in the present 
modification the orthonormal basis \k = \ ) m) of the discrete series state space is somewhat 
different from the usually used one. The modified orthonormal basis is given by 

\k = |, m) = \n), m = n + <7, (4.47) 

where \n) {n = 0, 1, . . . , oo) is the number-state basis, so that m = a, a + 1, . . . , oo. Then 
the completeness relation is 

oo 

£ |i,m)(i,m| = i. (4.48) 

m=cr 

In the modified HP realization ( |4.43j ) the index m goes from a and not from zero, as it is 
customary. Thus we get, for a > 0, a generalization of the usual SU(1,1) discrete series 
representation. 

By using the identity resolution ( 4.4C ), we construct the Hilbert space of entire functions 
f(z,a), which are analytic over the whole z plane. For a normalized state |/) of the form 
( CT ), we get 

f(z,a) = <jT (\z\,<T) + (z*,a\f) = ^C^/)—-— (4.49) 

n=0 [71 -t a). 

and this state can be represented in the \z,a) + basis: 

\f) = I d»(z,a) f {z *' a) \z,a) + . (4.50) 
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The orthonormal basis u n (z, a) in the Hilbert space of entire functions can be chosen corre- 
sponding to the number-state basis: 

~n+cr 

u n {z, a) = ^T (\z\,a) +(z*, a\n) = - (4.51) 

[Tl ~ r (J J • 

The generators K±{cr) and K?,((j) act on the Hilbert space of entire functions f(z,a) as 
linear operators: 

* + w=* a-m=^ + 4, *m=4A (452) 

We obtain, for example, 

K 3 (a)u n (z,a) = (n + o + ±)u n (z,cr), 

K + (a)u n (z,a) = (n + a + l)u n+ i(z, a), (4.53) 

K-(a)u n (z,a) = (n + a)u n ^i(z, a). 

We see that the unusual features of the modified HP realization manifest in the fact that 
the vacuum property, 

K-{a)uo(z,a)=Q, (4.54) 

must be included as an additional restriction. This fact is related to the use of the antinormal 
ordering, that is introduced as an additional restriction in order to restore the unitarity of 
the exponential phase operators. 

2. Statistical and phase properties 
The photon- number distribution of the \z,a) + states is 

1 \ z \2(n+a) 

Pi + \z,a) = \(n\z,a) + \ 2 = 1 (4.55) 

T {\z\,a) [{n + ay.} 2 

After some algebra we obtain 

1 o{\z\,cr) 

T {\z\,a) T (\z\,a) [T{a)\ 2 

where we have defined 

<t-2 |_|2m+l 



and the T (|,2|,cr) is defined by Eq. ( |4.38|) . The intensity correlation function is 
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9%{ 



a 2 + \z\ 2 + a- (2a + l)\z 



Tx{\z\,a) 



+ 



1 


2 


2<T 


To{\z 


TO] 2 



O" + \Z\ 



,T 2 (\z\,a) 



2a\z 



Ti(N,cr) 

T o(kl,^) 



(4.59) 



Only for a = and a = 1 the photon statistics is sub-Poissonian for all values of |z|, i.e., 
the gf\(+) is less than unity in the whole z plane. For a > 2, the g^l(+) is greater than 
unity while \z\ is less than a specific value depending on a. When \z\ exceeds this value, the 
intensity correlation function becomes less than unity. In the quantum limit \z\ <C 1, we get 



W£2 



<$(+) 



\z\ 



(o 2 -2)\z\ 



(a + 1)2 

12 



2 



(a + l) 4 (a + 2) 2 ' 
(3cr 2 + 4cr)|^| 4 



(a + 1) 2 ( CT + l) 4 (a + 2) 2 ' 

<7 + l 

a + 2 



(4.60) 
(4.61) 
(4.62) 



For a = 0, we return to the known value ~; for (7 = 1, the <?2 2 2(+) tends to |; for a > 2, the 
limiting values of the g z 2 l(+) are greater than unity. For very large values of a (a ^> 1), the 
intensity correlation function ( |4.62j ) tends to the thermal value 2. In the limit \z\ ^> a, we 
get 

1 



(n) 



(+) 



a 



I |2 , 2 



2o-(|z| 



1), 



1 - 



2\z\ 



(4.63) 
(4.64) 
(4.65) 



Thus, the \z, a) + states tend, in this limit, to have Poissonian photon statistics in the same 
way as the BG subcoherent states and the \z,a)- states. 
The phase representation function of the \z,a) + states is 



Q(+\z,a;9) 



Sad 



T o(M,o") 



(7-1 



exp (ze * e ) - X! 



ze 



m=0 



ml 



(4.66) 



When a is of order of unity, we can refer to the \z, a) + states as philophase states. However, 
when a increases, more and more additional terms are included in the ®( + '(z, a; 9) function, 
so that the \z,a) + states become unsuitable to be called philophase states. The phase 
distribution function is given by 
1 

~ 27T 1 



Q 



(+). 



Z,(7\ 



lew; 



cr-l 



ri+m 



2irT (\z\,<T 
- exp ^ze~ l9 ^ 



exp [2\z\ cos(9 — (p z )] + 

n,m=0 

cr-l 

- exp (z*e ie ) 



cr-l 



re 



nlml 



exp [i(n — m){9 — <p z )] 



ze 



-w 



m=0 



ml 



m=0 



ml 



(4.67) 
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The larger values of a are, the flatter the Q^(z,a;9) function is, and the worse the phase 
of the state is defined. By using the phase distribution ( (4.67] ), we calculate numerically the 
phase variance *(A0)^ CT (+)*. For \z\ — > 0, the *(A0)^ CT (+)* tends to the random value tt 2 /3. 
For 1 2 1 3> cr, the phase variance tends to zero. For given \z\, the larger values of a are, the 
more uncertain the phase of the state is. 

We study the number-phase uncertainty relation of the \z,a) + states by calculating 
numerically the V^ + \z, a) function, defined according to the general expression ( |2.72| ). When 
\z\ — > 0, the V^ + \z, cr) function tends, for any cr, to the standard random-phase value 7r 2 /12. 
While \z\ increases, the V (+) (z, a) function at first grows, reaches a maximum at a specific 
value of \z\ depending on a, and then decreases tending, for \z\ ^> a, to the limiting value 
|. For a = 0, the maximum is already at \z\ = 0. The larger values of a are, the higher 
the maximum of the V (+) (z,a) function is, and the slower an equality is achieved, as \z\ 
increases, in the number-phase uncertainty relation. 



3. Contraction to the Glauber CS 



Consider the boson annihilation operator a = e^yh acting on a photon state with mean 
photon number (n) and photon-number variance (An) 2 . Defining 5h = n — (n), we can 
write h = (h) + 5h. We choose a state for which An -C (n). Then we can approximate the 
annihilation operator a acting on such a state: 

5h 



a = e^J (n) + 5n ps e^J (n) 1 + — - = — -=e^{h + in)). (4.68) 
V V V 2 \ n )J 2 J in) 



We see that in the described case the annihilation operator a can be approximated, up to a 
numerical factor, by the operator Z(a) = e % ^(n + cr), 

1 Zip), o-=(h), An « (w). (4.69) 



2^ 



Therefore, the \z,a) + eigenstates of the operator Z{a) are, in the considered case, an ap- 
proximation of the Glauber CS which are the eigenstates of the a, 

\z,a)+ m \a), An « (n). (4.70) 

By comparing eigenvalues, we get 

a « -^=. (4.71) 

The mean photon number of the Glauber CS is (h) a = \a\ 2 . We suppose that in the 
considered case the mean photon number (h)yj of the \z,a) + states is approximately equal 
to the (h) a , and this gives a = {n)yj ~ \a\ 2 . Then we find from Eq. ( |4.71| ) 

\z\ w 2\a\ 2 « 2cr. (4.72) 
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The number variance of the Glauber CS is (An) 2 = |a| 2 . Then the condition An <C (n) is 
satisfied for |a| ^> 1. The conclusion is that the \z,a)+ states contract to the Glauber CS 
\a) for a ~ \z\/2 ~ \a\ 2 provided that \a\ ^> 1 (the classical limit). 

This result can be verified by calculating statistical properties of the \z, a) + states with 
\z\ = 2a. For large values of a, we would get, according to the discussed contraction, 

(n>W « (An)L(+) » a. (4.73) 



This is confirmed by numerical calculations. By taking \z\ = 2a in Eqs. ( 4.56|) and (|4.57 ), we 
find that the relations {h)^} /a and (An) 2 , a (+) / a quickly tend to unity, as a increases. The 
difference from unity for (n)y~J/a is about 5 x 1CT 3 , as a goes to 50, and for (An) 2 a (+)/a 
it is about 10~ 5 , as a goes to 20. Also, we calculate numerically the intensity correlation 
function g^l(+) of Eq. ( |4.59| ) for \z\ = 2a. This function quickly tends to the Poissonian 
value 1 of the Glauber CS. The difference of the gfl(+) from unity is about 5 x 10~ 4 , as a 
goes to 20, and it is about 10~ 4 , as a goes to 50. 



V. DISCUSSION AND CONCLUSIONS 

The SU(1,1) CS \k,C) and the BG subcoherent states \k,z) can be considered as two 
possible modifications of the familiar Glauber CS for the SU(1,1) Lie group. However, these 
two modifications lead us to the states with very different statistical and phase properties. 
The CS \k,C) have wholly super-Poissonian statistics, while all the BG states \k,z) are 
antibunched. With increase of the Bargmann index k, the photon-number distributions of 
the \k, C) and \k, z) states move from opposite sides to the Poissonian distribution. Another 
difference between the \k,() and \k,z) states is the behavior of statistical properties with 
change of the photon-excitation strength. As mean photon number increases, the relative 
photon-number uncertainty An/(h) tends to zero for the BG states, but it approaches a 
nonzero limit depending on k for the CS \k,C). Analogously, for large excitations, the 
intensity correlation function tends to unity for the BG states, but it does not depend 
on the excitation strength for the SU(1,1) CS. 

Phase properties of the \k,C) and \k,z) states show opposite behaviors with change of 
k. As k increases, the phase distribution Q(9) becomes narrower for the SU(1,1) CS, and it 
becomes flatter for the BG states. The behavior of the phase variance *(A0) 2 * with change 
of the photon-excitation strength is similar for all discussed types of states. In the quantum 
limit (near the vacuum) the phase variance tends to its random value 7r 2 /3, while in the 
classical limit (large excitations) it approaches zero, as must be for a classical wave with 
the perfectly defined phase. Phase properties of a photon state can be studied by using the 
cosine or sine variances as well as the variance of the Hermitian phase operator <p. Results 
are equivalent for using any phase-related observable, since optical phase is unique in the 
antinormal ordering formalism. 

An important problem discussed in the present paper is the uncertainty relation between 
the number and phase variables. Again, the choice of phase-related observable for studying 
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this uncertainty relation is arbitrary. We have seen that by using the phase variance *(A0) 2 * 
and the cosine variance *(AC) 2 * we obtain the same information. The properties of the 
number-phase uncertainty relation are essentially different for the \k,() and \k,z) states. 
Near the vacuum (\(\ — > or \z\ — > 0) the uncertainty function V tends to the random- 
phase value 7r 2 /12 for all types of states. However, in the limit of large excitations the 
behaviors of the SU(1,1) CS and the BG states are absolutely opposite. In this limit the 
V(fc, C) function of the CS \k,() blows up, according to the specific statistical properties of 
these states, while the V(k, z) function of the BG states \k, z) tends to its minimal possible 
value \. The larger values of k are, the slower the \k,z) states bring an equality to the 
number-phase uncertainty relation, as photon-excitation strength increases. In contrast to 
that, the CS \k, Q provide an approximate equality in the number-phase uncertainty relation 
for large values of k, in an intermediate range of excitations. With all that, the BG states 
\k,z) have a simple representation in the coherent-state basis \k,(), and vice versa. 

The case of the Bargmann index k equal to \ is interesting from two points of view. 
Firstly, this case is related to a simple type of intensity-dependent coupling in the Jaynes- 
Cummings model Hamiltonians. Secondly, the photon states associated with the HP SU(1,1) 
realization have, for k — |, special phase properties. Phase-state representation function 
Q{0), defined generally as Fourier series, can be converted into a relatively simple functional 
form for the philophase states ||,C) an d \\, z )- Moreover, the SU(1,1) CS ||,C) also are the 
eigenstates of the exponential phase operator e*^, just like the phase states \9). From the 
other hand, the phase properties of the BG states \\,z) can be generalized by introducing 
eigenstates \z, a) of the operator Z(a) = e"^(n + a). For integer a < 0, we find the class of 
generalized philophase states \z,a)-. All of them have the same phase distribution function 
as the BG states \\,z). The philophase states \z,cr)- are antibunched, and in the quantum 
limit \z\ — * they tend to the number states \n = \a\). In the case of integer a > 0, we find a 
modification of the HP SU(1,1) realization with k — \. Statistical and phase properties and 
the number-phase uncertainty relation of the states \z,a) + have interesting features. For a 
of order of unity, the \z,a) + states are close to be described as philophase states, while for 
(j>l, they contract to the Glauber CS \a), provided that \z\ ss 2a and a w z/(2y/a). 
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